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In discrete energy structure of nanoscale magnets, nonadiabatic transitions at avoided level 
crossings lead to fundamental processes of dynamics of magnetizations. The thermal environ- 
ment causes dissipative effects on these processes. In this paper we review the features of the 
nonadiabatic transition and the influence of the thermal environment. In particular we dis- 
cuss the temperature independent stepwise structure of magnetization at very low temperatures 
(deceptive nonadiabatic transition), the alternate enhancement of relaxation in the sequence of 
resonant tunneling points (parity effects), and processes caused by combinations of nonadiabatic 
transitions and disturbance due to external noises. 

KEYWORDS: Nanoscale magnets, nonadiabatic transitions, Landau-Zener-Stuckelberg transitions, dissipative envi- 
ronments 



§1. Introduction 

Hysteresis phenomena of ferromagnets have been one 
of the most interesting problems in the magnetism and 
statistical physics. Mechanism of the coercive force has 
been investigated by studying the processes that lead, tp 
the critical nucleation and motion of the domain wallBB 
From the point of view of free energy of the system, the 
hysteresis phenomena have been discussed in terms of 
the relaxation process of the metastable state to the true 
equilibrium state in the picture shown in Fig. 1 (a). 
Usually the end point of hysteresis is related to the so- 
called spinodal point where the metastability disappears 
(see Fig. 1(b)). 

However there is some probability of relaxation from 
the metastable state B to the equilibrium state A. For 
example at finite temperatures the probability of ther- 
mal excitation to the top of barrier C is proportional to 
exp(— AEs/hsT) and thus the relaxation rate pth of the 
metastable state through this activation process is given 

by 

1 -AEg 

Pth = — e fc B T , (1.1) 
To 

i.e., the Arrhenius law. At low temperatures the relax- 
ation time t = 1/pth of this process diverges exponen- 
tially. 

The shapes of the free energy in Fig. 1 are given by 
the mean field theory, which gives a good intuitive pic- 
ture of the metastability. Here it should be noted that 
AEb should be a microscopic quantity. Because in a 
bulk system AEb is of order the system size, the activa- 
tion rate vanishes, i.e., pth = 0. For a more quantitative 
understanding, we must look at the system microscopi- 
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Fig. 1. Potential pictures of relaxation: (a) metastable (b) spin- 
odal. 



cally. There the system is not uniform and we have to 
consider a microscopic break through of the metastable 
state. Such break through occurs as a process of creating 
the critical nucleus.ET For this microscopic process, Figs. 
l(a)-(b) represent effective potentials of the size of the 
nucleus. The relaxation time of metastable states has 
been classified according to the size of system and the 
generating rate of the critical nuclei. There are two re- 
gions, i.e., a single nucleation region (stochastic region) 
and a multi nucleation region (Avrami region) Jj' 

When the size of the magnets becomes smaller than 
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the width of the domain wall, the nucleus can not be 
defined. In such cases, the magnetization of the sys- 
tem changes uniformly and this process of breakdown of 
the metastable state is called "coherent process". Re- 
laxations in this situation also have been studied exten- 
sivelyD-* 

It has been pointed out that quantum fluctuation may 
play an important role in such small systems. In order 
to detect such quantum processes, several experiments 
have been proposed.™ However distribution of particle 
sizes prevents to analyze their processes in simple ways.™ 
Studies on single magnetic particles have been also per- 
formed but clear evidence of quantum processes has not 
yet been found .LP 

In this-xespect. nanoscale molecular magnets such as 
MnufjHBBBfci Fe 8 @>@@> and V 15 EP etc. are 
more promising. These molecules Mn 12 and Fes consist 
of small number of atoms. The low energy state of the 
system is represented by an effective S = 10 spin. Be- 
cause interactions among molecules are very small, each 
atom can be regarded as a S — 10 single spin. The 
Hamiltonian of the spin is generally given by 



H = - 

where S z = -10,-9, 



-DSi - HS Z 



(1.2) 



causes the quantum fluctuation, such as S x , S% — Sy 



10 and Q denotes a term which 

or 

(S + ) 4 + (>S~) 4 . In these systems the energy levels as a 
function of the field have a discrete structure (Fig. 2(a)). 
There we expect an explicit quantum mechanical dynam- 
ics. 

Due to the term Q, a small energy gap is formed at 
each crossing point as shown in Fig. 2. This structure 
is called avoided level crossing. When the field is swept 
through an avoided level crossing point, so called nona- 
diabatic transition occurs. Nonadiabatic transition plays 
important roles in microscopic quantum dynamics such 
as level dynamics of semiconductor, chemical reaction 
and optics. Nonadiabatic transitions in various cases 
have been reviewed by Nakamura.ES' 

In the uniaxial magnets, the nonadiabatic transition as 
shown in Fig. 2(b) occurs. Here the population coming 
in from the channel A is scattered to the channels B and 
C with probabilities p and 1 — p, respectively. Here the 
channel B is the ground state. Thus the scattering to B 
is an adiabatic change and corresponds to the tunneling. 
On the other hand, 1 — p is a probability to jump up to 
the channel C. This process corresponds to the case to 
stay in the metastable state (un-tunneling) . 

The states of the channels A and C are the same state 
when there is no quantum fluctuation i.e., Q = 0. This 
unperturbed state is called the diabatic state. These 
channels A and C have similar states even in the pres- 
ence QJLQ. The probability. p was studied by Landau,EJJ 
ZenerEf and StuckelbergEl (LZS) and is given by 

tt(AE) 2 



p = 1 — exp 



(1.3) 



2chg^Am j 

where AE is a gap at the avoided level crossing and 
Am is the difference of magnetization of the levels, c 
the speed of the sweeping field c = dH/dt. Thus the 
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Fig. 2. Energy structure of S = 10. Uniaxial magnet as functions 
of the external field H . (a) global structure and (b) avoided level 
crossing. 



tcri2T_au«AjTT.c is the changing rate of the Zeeman en- 
ergy .E3I3E5' The probability p plays an important role 
in quantum mechanical relaxation of the present system. 

In this LZS type nonadiabatic transition the transition 
occurs only in the vicinity of the crossing points. The 
first property is an essential ingredient for the relaxation 
at discrete points of the magnetic field in the molecu- 
lar magnets.tfP Making use of this dependence we can 
estimate AE from the change of magnetization AM . 
Let the crossing levels have magnetizations of the cor- 
responding adiabatic states m and m' . The change of 
magnetization is given by 



AM = pm! + (1 — p)m 



■ p(m' — m). (1-4) 



In the cases where AE is observed by other methods such 
as AC-susceptibility, this sweeping dependence of AM 
would give a method tOpCpnfirm that AE really comes 
from the tunneling gap.E9 However such confirmation 
has not yet been observed in small magnetic particles 
such as magnetic dots and ferritin. On the other hand 
the two characteristics have been observed in molecu- 
lar magnets at least qualitatively. As a more peculiar 
property of the nonadiabatic transition, it has also been 
pointed out that the application of an alternate field at 
the resonant points will cause a nontrivial oscillation of 
magnetization due to phase interference.^! 
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In order to study experimental data quantitatively, 
we need to incorporate the effects of environments. Ef- 
fects of noise oh. the LZS transition have been studied 
by Kayanuma.B3 The dynamics of molecular magnets 
in dissipative environments has also been investigated 
quite extensivelyrZP 

A thermal bath causes enhancement of the relaxation, 
e.g. the thermally assisted resonant tunneling, where 
resonant tunnelings of excited states play an important 
role. 

In this paper we study the characteristics of the reso- 
nant tunneling affected by thermal disturbance. In par- 
ticular, we study the effect of environment at very low 
temperatures such that relaxation process does not de- 
pend on the temperature. Even at such low tempera- 
tures, contact with the bath causes relaxation between 
levels which have magnetizations of the same sign. These 
states belong to the same valley in the potential picture 
of Fig. 1(a). The energy barrier does not exist between 
them. It is found that the relaxation between them easily 
occurs even with very weak disturbance. In pure quan- 
tum mechanical motion, transitions between these levels 
are almost prohibited except near the avoided level cross- 
ing points. Thus the magnetization curve with the dissi- 
pative effect is different from the one of a pure quantum 
case, although it does not depend on the temperature. 
We call such a process "deceptive nonadiabatic transi- 
tion" E3> 

At higher temperatures, excitation levels begin to con- 
tribute to the relaxation phenomena. At higher temper- 
atures, alternate enhancements of relaxation at resonant 
points are observedjf which is called 'the parity effect'. 
We consider the mechanism of such an alternation in a 
view point of nonadiabatic transition of excited states 
and hnd it as a general property of resonant tunnelingjQf 
excited states reflecting the structure of energy levels.EEP 
We also discuss the \/i-dependence of initial decay at 
resonant points. Furthermore we study various cases of 
the LZS process in fluctuating random environments. 

§2. Numerical Method 

The most standard method to study quantum dynam- 
ics in dissipative environments is the quantum master 
equation (QME) which describes the equation of motion 
of the reduced density matrix of the system p(t) which is 
derived by tracing out the degrees of the freedom of the 
environment from the density matrix of the total system. 
The total system consists of a system 7is, .a-thermal bath 
H.B and an interaction between them Tii-ES 

H = H S +H I +H B . (2.1) 

The reduced matrix is given by 

p (t) = Tr B e-^. (2.2) 

We have the following equation of motion for p(t) in 
the limit of weak coupling, assuming that the correla- 
tion time of the bath variable is very short (Markovian 
approximation) 

±p(t) = \[H,p(t)]+Tp(t), (2.3) 
at in 



where T is a linear operator acting on p{t) . This equation 
has been used to study quantum dynamics of optical 
process, etc. In most cases Tp has the so-called Lindblad 
formtiiP 

Tp = A ] Ap + pA^A + A^pA + Aprf , (2.4) 

where A is an operator of the system. However in mul- 
tileveled phenomena Tp has a more general form. 

In the cases where the bath consists of an infinite juim- 
ber of bosons, a general expression can be derived.^ 

^ = -i [H, p(t)} - A ([A, Rp(t)} + [A, Rp(ttf) , 

(2.5) 

where 

(k\R\rh) = a E ~ k ~ Efh )np(E- k - E*){k\X\m), 

C(w) = I(uj)- I(-uj), and np(u>) = — - ■ 

Here (3 is the inverse temperature of the reservoir 1/T, 
and we set Ti to be unity. \k) and |m) are the eigenstates 
of H with the eigenenergies Ej, and Em, respectively. 
I(ui) is the spectral density of the boson bath. Here we 
adopt the form I(ui) — Iqoj . When a = 1, it corre- 
sponds to the so called Ohmic bath and when a — 2, 
it corresponds to the phonon bath (super-Ohmic). As a 
more realistic bath for the experimental situation at very 
low temperature, we may take the dipole-field frorruithcr 
molecules and interactions with the nuclear spins. B A 
is an operator of system which is attached to bosons 
of the reservoir linearly, representing the interaction be- 
tween the system and the thermal bath. In the present 
study we take A = \ (S x + S z ) . The relaxation process 
depends on the form of A. Generally the coupling with 
the transverse component A = S x is more efficient than 
that of the longitudinal one X = S z for the relaxation. 
Detailed comparisons among choices of the form of the 
coupling will be presented elsewhere. 

For strong noise caused by fluctuating forces we can 
simulate quantum dynamics by solving the Schrodinger 
equation in random fields .ES 

§3. Quantum Dynamics in Dissipative Environ- 
ment 

3.1 Deceptive nonadiabatic transition 

In the lowest avoided level crossing point (— S, S) the 
change of magnetization AM is given by ( |1.4| ) . However 
at higher crossing points (m, m') with ml < S, the popu- 
lation scattered from m to ml is found to decay easily to 
the ground state, i.e., ml — > S, even when the dissipative 
effect is so small that the population at the metastable 
level of m hardly decays. This difference can be easily 
understood from the intuitive picture of Fig. 1(a). That 
is, the relaxation in the same valley, i.e., ml — > S, is easy 
while the relaxation over the barrier m — > S is hard . In 
this situation, we can not apply the relaxation (|l.4| ) di- 
rectly to estimate the LZS probability p. However we can 
still estimate p using AM because the relaxation from 
the level of m occurs with the LZS probability and the 
relaxation to the ground state occurs in a rather short 
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time. Taking these points into account, we modify the 
relation ( |l.4| ) by replacing the final magnetization rn! by 
S: 



AM = pS — (1 — p)m — m — p(S — m). 



(3.1) 



In order to confirm these processes we performed simula- 
tions using the QME. First we confirmed that relaxation 
from the metastable point is unlikely to occur when the 
coupling between the bath and system is weak and the 
temperature is low. On the other hand, a fast relax- 
ation is observed between levels with magnetizations of 
the same sign, which are in the same valley in Fig. 1(a). 
Furthermore when we sweep the field we find a step- wise 
magnetization curve whose step heights do not depend 
on the temperature but are definitely different from the 
pure quantum case. In Fig. 3, we show an example of 
magnetization process for T = 0.1, T = 0.5 with very 
small effects of environment (A = 0.00001: a solid line) 
and that of pure quantum system (A = 0: dashed line). 
Both of them are not temperature dependent within this 
temperature range. We call this stepwise structure in dis- 
sipative environments 'the deceptive nonadiabatic tran- 
sition'. We find that we can correctly estimate the pure 
quantum transition probabilities using the relation (3.1). 
Thus even at very low temperatures the effect of the en- 
vironment can not be excluded, but quantum mechanical 
processes and dissipative effects due to environments can 
be disentangled, and the information on the LZS proba- 
bilities can be extracted. 

For the phenomena described above, the existence of 
the environment is important but the detailed nature is 
not important as long as it leads to fast relaxation to the 
ground state. If the environment causes a change of LZS 
probability, which would be possible when the sweeping 
rate is very slow, further consideration is necessary.Er 




dent phenomena. These processes would depend on the 
detailed characteristics of the bath and the ways of the 
coupling between the system and the bath. Therefore 
general description is difficult. However here we point 
out a general property of relaxation under these condi- 
tions. 

As a characteristic of resonant tunneling at rather high 
temperatures, it has been observed that the amount of 
relaxation at the resonant points changes alternatively.EP 
Along a diabatic line, the energy gap increases monotoni- 
cally as the difference of magnetizations | m — m! | of levels 
decreases. Thus the transition probabilities at the reso- 
nant points increase monotonically. In a perturbational 
treatment the energy gap depends on the difference asEJ 



0.85 



Fig. 3. Magnetization processes in pure quantum dynamics (a 
dashed line) and a weak dissipative environment(a solid line). 



3.2 Parity effect 

At higher temperatures, excitation levels begin to con- 
tribute to the relaxation leading to temperature depen- 



AiJoc — 
D 



\m—m' | 



(3.2) 



Thus we have to consider a mechanism of the alternate 
enhancements. Here we interpret it from the view point 
of resonant tunneling of excited state. The transition 
probabilities at resonant points with the same value of 
\m — m'\ are nearly the same. Those points are located 
at the same horizontal level in Fig. 2. For example the 
values of p given by (1.3) for the case of T = 0.45 with the 
sweeping speed c = 0.0001 at the points, (—8, 5), (—9, 4), 
and (-10,3) are 0.91, 0.64 and 0.99, respectively. On 
the other hand, those at (—8, 6), (—9, 5) and (—10, 4) are 
0.72, 0.037 and 0.01, which are very small. Thus most of 
the population at the levels m = —8, —9 and —10 decays 
at the former points. These decays cause enhancements 
of relaxation at H = 0.3,0.5 and 0.7, which gives the 
parity effect. In Fig. 4, we show the magnetization of 
this case with its time derivative. 

Because the energy structure shown in Fig. 2 is gen- 
eral for uniaxial magnets, we expect that the alternate 
enhancement of relaxation, i.e., the parity effect, is a gen- 
eral property of resonant tunneling in the thermal envi- 
ronment. We have also pointed out that if we change 
the sweeping rate the enhanced sequence is shifted. For 
example if we sweep much slower, the probabilities at 
(—8, 6), (—9, 5) and (—10,4) become large and popula- 
tions on the lines decay there, which causes the shift of 
the enhanced sequence at H = 0.2,0.4 and 0.6. 



3.3 N on- exponential decay at the resonant point 

The magnetization which is initially polarized upward 
decays rather fast at a resonant point. Here the field is 
set at this point and is not swept. There are several paths 
for the magnetization to relax at this point. First let us 
consider relaxation by the nonadiabatic transitions at 
the lowest resonant point. Because at the resonant point 
the energy gap is very narrow, the field fluctuates around 
the point as shown in Fig. 5(a). If we regard the motion 
of the field as a Brownian motion, it is known that the 
number of the times the field crosses the resonant point 
is proportional to ^/i, i.e., thc-recurrence time of one- 
dimensional Brownian motionBa) 

At each crossing, the population moves to the other 
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Fig. 4. The parity effect: Magnetization process (a solid line) and 
its time derivative (a dashed line). 
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Fig. 5. (a) Random field near the resonant point and (b) Relax- 
ation process through the excited state. 



branch by the LZS transition probability 

Pi = 1 - cxp( ) ~ — — . (3.3) 



4c,; 



In a pure quantum mechanical process, quantum me- 
chanical interference occurs among the transitions SB 
But assuming a fast decoherence, the total transition 
probability is expected to be given by 

Ptotai(*) ^ J2 Pi ~ ^(~) = aVi. (3.4) 

Thus we naturally expect that the magnetization decays 
as 



AM = M (l - 2aVt) 



(3.5) 



at the initial stage. 

For a longer time scale, the field does not fluctuate 
freely but is confined near the resonant point. Thus 
Ptotai for long time is proportional to t. Therefore the 
magnetization decays in exponentially 

PtotalM 



This mechanism described above may give the simplest 
explanation of the y/i behavior. A-pore detailed analysis 
has been given in the reference,0 taking into account 
the explicit nature of the fluctuation of the external field. 

3-4 Modification of the transition rate 

When the amplitude of the external disturbance is 
strong, we have to use another estimation of the tran- 
sition probability. Kayamura and Nakayama have inves- 
tigated LZS transition in fluctuating ficld-and obtained 
expression of the transition probability.© In the case 
where the sweeping rate is slow and the transient time 
through the resonant point is much longer than the phase 
coherence time, then they obtained the transition prob- 
ability as 

P = Psd = \{l - exp(~TT(AE) 2 /2c%^ B Am)). (3.7) 

On the other hand, when the transient time is very short, 
transition probability does not change from that of the 
pure LZS transition 

P = Plzs = 1 - exp(-7r(A£:) 2 /4c%^BAm). (3.8) 

They confirmed such dependences by numerical simula- 
tion. 

Next, let us consider a case where excited levels con- 
tribute to the relaxation. If the frequency of contacts 
between the system and the bath is high, a tunneling 
through the excited state would enhance the relaxation 
rate-even the population at the excited state is very lit- 
tle£3> 

Let us consider the case where the LZS transition prob- 
ability at the lowest level po is very small and the one at 
the excited level p\ is of the order 1. Let us consider a 
case where a state is excited to the exited level. In the 
off-resonant region this excited state decays to the orig- 
inal state very rapidly. The population at the excited 
state tie is determined by the balance equation 



(3.9) 



where uq is the population in the ground state and 
Re^g an d i?G_>E are transition rate from the excited 
state to the ground state and vise versa. Although 
riE = Rg^e/Re^git-g is very small, Re^g an d i?G— »e 
themselves can be very large. 

At the resonant point, small fluctuations of the field 
would cause the crossing point as we see in Fig. 5(a). 
Thus the population pumped to the excited state A' can 
be transferred to the state B' and then it decays to B 
instead of A. (Fig.5(b)) This path (A-> A'-> B'-> B) 
becomes dominant when p\ po and fluctuation of the 
field is rather fast. The opposite path (B— > B'— > A'—* A) 
is also larger. In the initial stage most of the population 
is at A and therefore the population moves along the 
former path. Thus the effective transition rate A — > B is 
enhanced very much. 



(3.10) 



M(t) 



-t/r 



t 



(3.6) 



where v is the frequency of the crossing. For a short 
time v is proportional to \ft and for a long time it is 
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proportional to t as we saw above. 

To study the transition probability in this case, the 
correct information for Rq— >e and v is necessary, but it 
is generally difficult. However, if we could estimate these 
quantities from the enhancement, it would yield detailed 
knowledge of the bath. In nanoscale molecular magnets, 
it would be possible to study such a detailed property, 
which is a very interesting research area in the future. 

§4. Summary and Discussion 

Nanoscale molecular magnets display several phenom- 
ena which originate from explicit quantum mechanical 
transitions between discrete levels. In this paper we 
studied effects of dissipative environment which smears 
out the pure quantum processes. So far the relaxation 
processes have been a kind of black-box and have been 
treated only phcnomenologically. But it would be pos- 
sible to begin to study explicit processes of relation in 
nanoscale magnets because of their simple form. 

So far we studied the S = 10 spin representing the 
low energy structure of magnetic energy levels. Let us 
consider the structure of the full energy level. The molec- 
ular magnets have complicated structures. For example 
Mni2 includes 12 Mn molecules with many other atoms 
which have nuclear spins. Thus the dimension of total 
Hamiltonian is 5 8 • 4 4 x /, where / comes from the degree 
of freedom of nuclear spins. This degree of nuclear spins 
causes random effects on each Mn atom. It would be an 
interesting problem to study how this random field on 
individual atoms causes changes of the energy levels at 
low temperatures. 

Even without the effects of nuclear spins, there are 
dipole-dipole couplings among the molecules which cause 
random noise on the whole molecule. Thus it should be 
taken into account even in the view point of S = 10 spin. 
Effects of this field are studiecLas 'feedback-effect' on the 
LZS process of magnetization^ It would be interesting 
to study natures of noises explicitly in nanoscale molec- 
ular magnets. 
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